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Abstract 

We motivate metrology schemes based on topological singularities as a way to build robustness 
against deformations of the system. In particular, we relate reference settings of metrological 
systems to topological singularities in the measurement outputs. As examples we discuss optical 
nano-position sensing (i) using a balanced photodetector and a quadrant photodetector, and (ii) a 
more general image based scheme. In both cases the reference setting is a scatterer position that 
corresponds to a topological singularity in an output space constructed from the scattered field 
intensity distributions. 
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I. INTRODUCTION 


Measurements are the way to access physical properties of a system and are indispensable 
in natural sciences. However, the design and assessment of measurement strategies, which 
is the primary goal of the field of metrology, can be a complicated task e.g., due to the 
presence of noise and the intrinsic complexity of the system at hand. The precision of 
a given metrology scheme can be assessed using metrics such as the Fisher information 
and the Cramer-Rao bound [I]. The Fisher information determines the maximum amount 
of information about the parameter that can be extracted from the system using a given 
measurement. The Cramer-Rao bound is an upper bound on the precision with which the 
parameter can be estimated and depends on the Fisher information. However, in practical 
experiments one is also interested in the robustness of the given scheme to the presence of 
noise and in its versatility , namely, its application to a broader set of systems. 

In this paper we focus on the robustness of metrology schemes, and our analysis also 
naturally leads to a general metrology framework for designing versatile schemes. Inspired 
by the use of topological methods to classify physical phenomena both in the classical and 
quantum regimes [2], to enable fault tolerant quantum computing [3], and for feature de¬ 
tection in image processing HHE!, our strategy is based on identifying topological features 
from measurement outcomes. We exploit the inherent stability and robustness of these 
topological features |7J, which make the metrology scheme resilient to a large class of de¬ 
formations, including stochastic noise and imperfections of the experimental method. We 
illustrate the approach in the context of nano-position sensing where we relate the reference 
position of the system to a topological singularity apparent in the space constructed from 
measurement outcomes. We will present some examples where we demonstrate that the 
topological singularity used as a reference position of the system is only revealed through 
the measurement and data analysis, i.e. it is not physically attached to the system. We 
also show that designing a metrology scheme based on topological features makes it robust 
to certain deformations of the physical system , by showing how the topological singularity 
survives imperfections of the elements used in the positioning experiment or deformations of 
the target object. In general, we claim that designing a metrology scheme on topologically 
distinct features of the system leads to versatility and robustness of the method. 

The paper is organized as follows. In Sec. [TT] we lay out the mathematical formalism 
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and explain the general idea behind our topological approach. Section |TTT] describes simple 
illustrations of our approach in the context of position sensing. We explain the familiar 
example of position sensing using a balanced photodetector, a quadrant photodetector, and 
then discuss a generalization of the image analysis at the output. Section [TTT] may be read 
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independently of Sec. ITT} We conclude with a discussion in Sec. IV 


FIG. 1. An abstract model of a metrology scheme consists of three topological spaces (i) 1Z the space 
of physical parameters to be measured which describe the system, (ii) A4 the space of the outcomes 
of the measurements of the physical parameters, and (iii) O the “order parameter” space. The 
spaces are equipped with continuous maps / : 1Z —>• M. and g : M. —>• O. The desired metrological 
information, e.g. the reference position in a position sensing experiment, is extracted from the 
order parameter. A topological metrology scheme is one in which the metrological information is 
extracted from a topological feature , such as a topological singularity, of the map g o f. 


II. TOPOLOGICAL METROLOGY FOR REFERENCE DETERMINATION 

A general metrological setup typically consists of three components which can be modelled 
as topological spaces as shown in Fig. [lj (i) 7 Z the space of parameters of the physical system, 
which are to be measured, (ii) M. the space of measurement outcomes corresponding to 
those parameters, and (iii) O the “order parameter” space—the image of a map g : M. —» O 
that translates the bare measurement outcomes to a quantity that provides the desired 
metrological information. We also denote by / : 7Z Ai an abstract continuous function 
that maps, in a given trial of the metrology, the physical parameters to a measurement 
outcome. 

In general, 1Z could be the configuration space of positions and velocities for a classical 
dynamical system, or a Hilbert space describing a quantum mechanical system. To be more 
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concrete here, we consider a generalized optical position sensing setup and in the following 
sections we will apply them to position sensing experiments. In this case 71 = the 
space of D —dimensional position vectors with coordinates (aq, x 2 ,..., xp) describing the 
location of the object. M. is the space of real scalar fields I(yi, 1/2, • • • , Vd'\ %i, ^2, • • •, xd) 
on M D \ and these fields typically correspond to the measured intensity distributions. No¬ 
tice that the dimensions (D 1 ) of the space where the intensity measurements are taken do 
not need to be the same as the dimensions (D) of the space of physical positions to be 
determined. In practice, M. D ' is discretized to the set of points {ii, • • • where each 

index if. can take values from if. = 1 ,,nk, with rq, = N, where N would represent 
the total number of image pixels, such that the continuous function is mapped to a dis¬ 
crete set: I(yi,y 2 ,---,yD r ,xi,x 2 ,...,x D ) -A {Ii u i 2 ,...i D ,{ x i,X 2 , ■ ■ ■ ,xd)}- We will focus on 
the particular goal of determining a reference position in the experiment, which may be 
required for example to align different components of the system or to reliably return to 
an initial position of the system. We require that the method should still yield a reference 
position in the presence of noise, such as intensity fluctuations of the input optical field, or 
electronic noise in the detectors, and possibly other deformations in the experiment, such 
as differently shaped target objects or differences in the efficiency and gain of the detectors. 
Here we advocate to introduce an order parameter space and the associated map g such 
that the reference position corresponds to a topological singularity in the composition g o f. 
To concretely demonstrate the topological analysis, we choose our order parameter space O 
to be the space of real unit vectors v = u(aq, x 2 ,..., xd) in the D dimensional Euclidean 
space E D with components denoted Vj. All vectors at non-singular points are normalised 
as JT v'j = 1. From now on, summation over repeated indices is assumed. We have chosen 
that the number of components of the vector v is the same as the number of dimensions in 
our reference space, D, because in that case the only possible stable topological singulari¬ 
ties that can appear will have dimension zero (points) [7j, which matches our aim to find a 
reference position. Let V be some region in K D and S its closed boundary: S = dV. At the 
non-singular points of v, define the differential D — 1 form 

uj = Vj 1 dvj 2 A dvj 3 A • • • A dvj D e^ 2 '"^ D (1) 

where e is the rank D totally anti-symmetric tensor. Consider the following quantity 

T(D/2) r 

" = i 
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with T(x) being the gamma function. Appendix [A] contains a formal derivation of the 
following important facts: The quantity /i can only have an integer value and is a topological 
invariant of the system. The actual value of /i can be calculated by summing the topological 
indices of singular points, p , of the system, by deforming the contour of integration, i.e. 
// = Up- Here p p can be formally defined as 

f ‘ r=I 2^£ weZ —{°t (2) 

where <f s surrounds a neighbourhood V p small enough that it only contains singularity p. 
The following are some illustrative examples of the use of this formalism. 

1 . D = 1 . A system whose position is mapped to an intensity field /^(aq), with A = 
{1,2}. Here, we have chosen D' — 1 and discretized this space so that N = 2. 
An order parameter which is sensitive to roots in the intensity function is u(oq) = 
g ° f( x l) = sgn(/ 2 (oq) — Ii(x\))y\. Over the closed interval V = [a, b] e M with 
boundary S = {a} U {6}, the index is 

/' = 9°f(b ) ~9° /(«))• 

This value will not change under deformations of the interval S provided a zero of the 
intensity is not crossed. 

2. D = 2 and D' = 0, so that N = 1. A system whose position is mapped to an intensity 

field I(xi, x 2 ). An order parameter which is sensitive to critical points in the intensity 
field is the normalised gradient v(x 1 ,x 2 ) = go f(x u x 2 ) = (l9 . J(a;i ' The 

invariant 

h = ^ § s (vid v 2 - v 2 dv i) 

= A/ dij ( dM. frWrfann 

2t tJS (djldjiy/ 2 \idjldjl) 1/2 (djldjl) 3 / 2 ) 

&2j ( ddil _ diIJPjIddjiy \ 

(djldjiy/ 2 \(djldjl) 1 / 2 (. djldjl ) 3 / 2 ) 

A d diIdd2tzhMdlL 
2tt JS {djldjl) 

Here the index is none other than the winding number. For example, if I{xi,x 2 ) = 
ax\ — bx 2 for a,b > 0 then choosing S to be the unit circle centred at the origin, 
/%,=( o,o) = — 1- 


5 













Having established the general criteria for a topological invariant determined by the unit 
vector held, we now cast the problem of measuring physical observables in that setting. The 
key observation is that we can perform measurements on a physical system and construct 
a unit vector held from those measurements in a way that singular points provide useful 
information, such as the centre of mass position or spatial symmetry of an object. In many 
circumstances, direct detection of a singular point may be difficult due to hnite resolution 
detectors, noise in the measurement, etc. However, indirect detection is always possible 
by computing the index // p over a closed path surrounding a region in which the point is 
proposed to exist, by Eq. (§ 

We wish to characterise the degree of stability of the index // p to noise and hnite resolution 
detection. As stated, the composition g o f dehnes a mapping from 7Z —» O. The quantity 

is invariant to continuous changes in / and g which deform the region V p but do not 
introduce new singularities that cross S p . Also, /j p is invariant under certain rescalings. For 
example, write the mapping g to the order parameter as a composition g = r o g where 
the image of g is the unnormalised vector held: v{x\,... ,xa) = g(I), and v = r(v). Then 
arbitrary rescalings v 3 —» XjVj, where X 3 > 0, will leave c o (Eq. 0 ) and hence //, p (Eq. (|2])) 
invariant. In fact if one only cares to detect the presence of a singularity, then the sign of /i p 
is irrelevant and it suffices to have X 3 all strictly positive or all strictly negative everywhere. 


III. APPLICATION TO POSITION SENSING 

A. Physical setup 

We now describe position sensing as a concrete example for a metrology scheme where the 
reference position is determined from a topological singularity. The physical system consists 
of an optical illumination held, a scatterer (e.g., a nanoparticle or a nanohole) that can be 
translated in the transverse plane, and a detection plane, as illustrated in Fig. [2] (a). The 
scatterer position is a vector in R 2 , which corresponds to the space 7 Z of the previous section. 
The optical intensity distribution is measured at the detection plane. Then, the space M. 
corresponds to all possible scalar intensity distributions in M 2 . Given that the scatterer 
can be moved at will, the aim is to reliably identify a reference position that we can easily 
return to. This will allow us to measure displacements from that position or simply use it 
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FIG. 2. (a) Physical system for optical position sensing: An incident light field scatters off the 
object in question, undergoes optional further transformations, and is measured in a plane behind 
the scatterer. The aim is to identify one specific position of the scatterer, which can be translated 
transversely, based on the measured intensity distribution, (b) Mapping from the M 2 space of 
scatterer positions (reference space TV) to unit vectors (in the order parameter space O) through 
the image analysis of the intensity distribution in space A4. The intensity image corresponds to a 
single scatterer position (xi,X 2 ), and thus we have dropped the dependence on these variables for 
the intensity, i.e. I(y u y 2 ] x lt x 2 ) = I(yi,V2) 

to stabilize the position of the scatterer. We then seek a reference position whose existence 
is robust to noise in the measurement. As we are in the linear optics regime, a small change 
in the illuminating field or in the shape or properties of the scatterer will result in a small 
change in the output field. This means that if we can find an idealized situation where 
the metrological scheme presents a topological singularity in the order parameter space, 
the structural stability theorems in the previous section will ensure that the topological 
singularity is also present when there are small deviations from the ideal situation. 


B. Example in 1 dimension: Balanced photodetector 

A common position sensing device used to measure the position of a scatterer along one 
axis is the balanced photodetector, illustrated in Fig. [3| One property that makes this 
method particularly popular is that it is quite robust under the effect of noise and other 
imperfections of the system. Such robustness was already apparent in the predecessor of 
position sensing devices [8]. We will now show that this method can be analyzed with the 
tools explained in the previous section. This will unveil that the main properties of this 
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FIG. 3. The balanced photodetector, a 1 dimensional position sensing device. The power incident 
on each of the two detector regions is amplified by a gain, and then the difference between the 
signals is evaluated. 

metrological method are founded on the identification of a topological singularity of the 
vector field constructed from the measurements. 

A balanced photodetector consists of two adjacent photodiodes with gain amplification 
factors Gl and Gr. For each scatterer position, we take a measurement of the power 
incident on the two regions, Pi, i G {L, R} which is amplified to produce output signals 
Si = P % x G{. A differential signal, AS = Sr — ,Sh, is calculated from each measurement. 
The ^-coordinate of the scatterer at which the signals are balanced, that is AS = 0, can be 
used as a reference. In order to formally follow the analysis in Sec. [TT| we will first consider 
the differential signal as a vector in E: V — ASy\. Then, we build the normalized vectors 
(v = R/jVj = sgn(AS')yi) in order to carry out the analysis in the order parameter. Now, 
the equivalence to the 1 dimensional case discussed earlier becomes evident. The map go f 
simply maps the x | -coordinate of the scatterer to a unit vector in E and the position where 
the scatterer balances the detector (AS = 0), appears as a topological singularity in the 
order parameter. 

Figure [4] illustrates the effect of deformations of the system on the measurement results. 
Shown are the only component of the unnormalised V and normalised v as a function of 
scatterer aq-coordinate for three different cases: (a) a case where the gains of the two regions 
are equal and no noise is present, (b) a case where the gains are unequal and still no noise 
is present, and (c) a case with equal gains but in the presence of stochastic noise with a 
bounded amplitude. Here we have considered the power gains as a perturbation of the 
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(a) (b) (c) 

G|_ — Gr G|_ < Gr G l ,G r noisy 



-1 - 0.5 0 0.5 1 -1 - 0.5 0 0.5 1 -1 - 0.5 0 0.5 1 


Scatterer x r coordinate Scatterer x 1 -coordinate Scatterer x,-coordinate 

FIG. 4. Typical output values V yi and v yi from a balanced photodetector with respect to the 
position of the scatterer in the aq-coordinate for (a) an ideal physical system, (b) unequal gains, 
and (c) equal gains but with added noise. 


physical system. Another component that could change is the probe light field. Since case 
(a) corresponds to an ideal physical system, its result can be regarded as the unperturbed 
reference position. The imbalance of the gains in (b) results in a perturbation that maintains 
the monotonicity of V yi but shifts the location of the single zero crossing—the singularity 
in v yi . The noise in (c) breaks the monotonicity of V yi and introduces further crossings in a 
way that the topological index fi remains invariant, which in the 1 dimensional case can be 
-1,0, or 1, depending on the sign of v yi at the boundaries of the interval in question (here 
fi — 1). Since the maximum amplitude of the noise is limited, the new crossings in V yi 
and hence singularities in v yi , lie within a bounded distance from the unperturbed reference 
location (corresponding to case (a)). In this case, it is no longer possible to pinpoint a 
single singularity as the reference, but we can still identify an interval that must contain the 
reference. This illustrates that while noise can impact the precision of the measurement, 
the existence of the reference position can nonetheless remain robust. This simple example 
should make our motivation for relating the reference position to a topological singularity 
apparent. 
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C. Example in 2 dimensions: Quadrant photo detector 


The quadrant photo detector (QPD) provides a simple extension of the position sens¬ 
ing scheme just described to 2 dimensions. A QPD consists of four detector regions with 
amplification gains Gi, i G {NE, NW, SW, SE}, and arranged as the four quadrants of 
the measurement plane. We perform a measurement of the power incident on the quad¬ 
rants, Pi which is amplified to produce output signals: Si = P, x G t . Two differential 
signals are calculated from each measurement, one horizontal difference of the signals, 
ASh = Sne + 'S'se — 'S'nw — 'S'sw, and one vertical difference, A S v = 5ne + <Snw — <Sse — Ssw- 
The scatterer position at which the signals from the four quadrants are balanced can be 
used as a reference point, which is equivalent to the condition AS'h = A S v = 0. This case is 
similar to the 2 dimensional case discussed in Sec. [TT] after we identify ASy, and A S v as com¬ 
ponents of a vector V = AS^yi + AS v y 2 in R 2 , which is normalised to obtain v = V/Vv ■ V. 
That is, the condition Ahj, = A S v = 0 corresponds to a topological singularity in the vector 
field v. In the idealised case where the gain factors of the quadrants are equal and both 
the scatterer and probe light have discrete four fold rotational symmetry (symmetry under 
rotations by it/2 about an axis through the origin perpendicular to the scatterer plane), it 
is clear that both differential signals are zero when the scatterer is placed at the symmetric 
position. On the other hand, this position corresponds to a singularity of the order param¬ 
eter, as argued in Sec. [TTJ In fact, such a reference position exists even when the light field 
and scatterer become asymmetric, and when the gain factors become unequal, since they 
correspond to continuous deformations on the order parameter space. This effect is equiva¬ 
lent to the one shown for the one dimensional case. The topological singularity is displaced 
due to the deformations in the gain or probe beam, but its structural stability avoids its 
destruction, unless the parameters change so radically that the symmetry of the system is 
completely lost. Again, a proper analysis shows that the QPD being robust, versatile, and 
easy to operate, which contributes to its popularity, stems from the topological features of 
the method. 
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D. Image-based nanopositioning using a topological singularity 

Finally, we describe a new, more general image-based position sensing approach, within 
the same physical setup described above, where the reference position is once again obtained 
from a topological singularity in the order parameter (space O) that is constructed from the 
output images. In a measured intensity distribution we define two points p[ and p 2 and 
consider the difference vector, V = P 2 — Pi, normalized to a unit vector v = V /\/v -V in 
M 2 . In this way, we map each position of the scatterer to a unit vector v in M 2 , see Fig. [ 2 ] 
(b). The vector field v corresponds to our order parameter as a function of the configuration 
space position. Our goal is to obtain a vector field v that changes continuously with the 
scatterer position, and that exhibits a topological singularity which can be used as a reference 
position. 

How can we choose the two points pi and p 2 to generate such a vector field? We can always 
keep one of the points, say pi, fixed. Then as long as p 2 shows a continuous dependence on 
the scatterer position, the vector field also varies continuously with the scatterer position. A 
quick analysis will show that this is a straight generalization from the BPD or QPD systems 
to an arbitrary large number of detectors, or pixels of a camera. Another, more interesting, 
option is to have both points depend (differently) on the intensity distribution. Care has to 
be taken when choosing the points to fulfil the conditions expressed in Sec. [TT| For instance, 
one could consider using the global minimum of the intensity distribution as one of the 
points. However, even when the intensity varies continuously with the scatterer position, 
the global minimum may undergo a discrete jump (e.g. the global minimum could jump from 
one trough in the intensity distribution to another distant trough that has approximately 
the same depth, even when the scatterer position is changed slightly). Such a choice of image 
analysis does not guarantee the desired property that the vector field changes continuously 
under a continuous change of the parameters. 

Here we consider defining the points Pi (i — 1,2) as weighted centroids, which satisfy the 
stated continuity requirement: 



( 3 ) 


where hi(I(yi, y 2 \ aq, x 2 )) is any continuous function of the intensity value at the coordinate 
(yi, y 2 ) in the image, when the scatterer position is (x 1 ,x 2 ). If the image is discretized due 
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( 4 ) 


to uniform rectangular pixelation, Eq. (|3j) is replaced with 

Pi = ' 52 h{Pj)(yi j yi + mjm), 

3 

where j labels the pixels, Pj = f A I(yi, 1 / 2 ', x 2 )dy 1 dy 2 with A, being the area of pixel j, 
and yip y 2 j are the mean y\,y 2 positions of pixel j. 

When considering a cylindrically symmetric light field and scatterer, the intensity dis¬ 
tribution has a rotational symmetry. As shown in the appendix, when the position of the 
scatterer (aq, x 2 ) is such that it lies on the symmetry axis, the cylindrical symmetry of the 
system enforces the weighted centroids in Eq. (|3]) to be zero, and the vector v in the order 
parameter has a discontinuity, Pi = p 2 = 0, provided a coordinate change in the image space 
to have the origin coincide with the symmetry axis. Thus, one can conveniently find a vector 
discontinuity in O in the presence of rotational symmetry. But what guarantees in this case 
that this vector discontinuity is indeed a topological singularity? By symmetry, a transla¬ 
tion of the scatterer between two positions equally distant from the center is equivalent to 
a rotation of the complete physical system, as illustrated in Fig. [5} This in turn leads to 
a rotation of the intensity distribution and of the vectors pi and p 2 . On the other hand, 
when the scatterer is placed at the central location, the two points coincide, resulting in a 
zero vector. Now, to prove that this scatterer position gives rise to a point singularity in the 
order space, choose a new scatterer position x that lies in the neighbourhood of this central 
position. Due to the symmetry breaking, the vectors pi and p 2 will no longer coincide, and 
the order parameter vector v will be well defined. As the scatterer traverses a concentric 
circle, v undergoes one complete rotation through 2n. Thus, the central point is a singularity 
with winding number p = 1. 

In a real experiment, perturbations generally break the cylindrical symmetry of the sys¬ 
tem. However, the same arguments used in the previous cases hold, and the topological 
singularity is robust under noise or perturbations of the system. Then, in the experiment, 
the topological singularity, which is related to the reference position, still exists as long as the 
corresponding perturbations in the vector field are continuous and sufficiently small. This, 
in fact, is the main motivation for relating the reference position to a topological feature 
obtained from the measurement outcomes. 

An interesting feature of this choice of image analysis is that it renders this positioning 
technology self-referenced. The reference position does not depend on the external devices 
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FIG. 5. Relationship between two optical systems where the scatterer has been translated. Left 
column: For a general scatterer and incident field, the translation of the scatterer can be replaced 
by a rotation of both incident field and scatterer, followed by a rotation of the whole physical 
system in the opposite direction. Right column: When the scatterer and incident field are both 
cylindrically symmetric and translation occurs between two positions of equal distance from the 
center of the incident beam, the same relationship applies. However, the rotations of the scatterer 
and incident field about their respective axes have no effect. Subsequent rotation of the complete 
system results in rotation of the scattered field pattern. 


used to measure it, as was the case with the QPD and BPD. In this case, the reference 
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position only depends on the incident beam and the scatterer. In the perfectly cylindrically 
symmetric situation it will coincide with the axis of symmetry of the incident beam and 
the scatterer. In real experimental situations it will correspond to the position of the stable 
topological singularity of the system. 


IV. DISCUSSION 

In this work we have presented a new way of analyzing metrological systems based on the 
geometrical properties of the experimental set-up and posterior analysis. These geometrical 
properties can be unveiled with the use of mathematical techniques in the field of topology. 
Then, the distinct features of the metrological system can appear as topological singularities 
in an order parameter derived from the measurement. We have shown that when such an 
approach is possible, the topological singularity serves as a distinct point in the system, 
which is robust under continuous perturbations. This can be used to set a reference point. 

Topological singularities present a structural stability under continuous changes of param¬ 
eters. Consequently, the existence of a reference point finked to that topological singularity 
will be particularly resilient to both external noise and the influence of small changes of 
parameters of the metrological set-up. We have exemplified all these features in the con¬ 
text of position metrology. Using the proposed approach, we have analyzed two common 
metrological systems used to provide a reference point (the balanced photodetector and the 
quadrant photodetector), and showed that most of the interesting features of the methods 
(including robustness and versatility) are connected to the fact that they are referenced to a 
topological singularity of the vector field obtained through measurement and data analysis. 
We have also proposed a new topological metrology method for position sensing, which offers 
an additional advantage that it does not rely on an external reference point. 

Topological concepts have already been used in image processing algorithms An 

important distinction is that the topological singularities in those references are at a different 
level compared to what was discussed in this paper. In the references, the vector field 
is obtained from a single image, and a topological singularity is identified in the image 
plane (the equivalent of (y 1,1/2) in Fig. [ 2 ] (a)). In contrast, we are taking an image for 
each coordinate in the reference space, which is the scatterer position (aq,^)- We do 
not seek to identify a point in the image plane as the reference point; each image only 
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provides one vector of the vector field. The topological singularity then resides in the vector 
field that is constructed in the analysis using multiple images. In this way, a point in the 
configuration space, i.e., a particular reference position of the scatterer, can be identified 
as the reference. We make this distinction by writing that our topological singularities 
are in an order parameter, constructed from measurement outcomes, as a function of the 
configuration space setting. We would like to stress that this concept is not limited to image- 
based measurements. The formalism presented is very general and could accommodate other 
metrological methods, such as measurements of magnetic fields, etc. The key ingredients 
are that the measurements have to be relative measurements, as most measurements are, 
and thus a reference point must be found, and that such reference point can be linked 
to a geometrically distinct point of the metrological system, in a way that when analyzed 
appropriately, it can give rise to a topological singularity in the order parameter space. 

As we have shown, the practical use of topological singularities in a vector field con¬ 
structed from measurement outcomes is key to understanding already existing metrological 
systems. The identification and analysis of topological features allows us to understand what 
makes such techniques and devices robust. Moreover, the analysis could also be applied to 
the design of new metrology schemes. An important point for further study is that from 
the analysis performed in this work, there is no immediate connection to the measurement 
precision. Consequently, optimal precision is not guaranteed by designing a technique in 
the way we have put forward. However, assuming that the precision of the measurement is 
limited by the precision in the location of the reference point, one could in principle use the 
theorems on the structural stability of the topological dislocations in order to give upper 
bounds to the precision of a topological metrology system for a given noise spectrum. On the 
other hand, we have discussed the robustness of the system. This is of practical importance 
because it affects the versatility of the method: Methods that use a topological singularity 
of an order parameter constructed from measurement outcomes as a reference are expected 
to be robust to deformations, and therefore, applicable to a number of physical systems. 


V. ACKNOWLEDGEMENTS 

This work was funded by the Center of Excellence in Engineered Quantum Systems 
(EQuS). G.M.-T is also funded by the Future Fellowship FF110100924. 


15 



Appendix A: Topological invariants 


The contents of this appendix can be found in most elementary textbooks in topology 
and other specialized texts m- We repeat it here for the convenience of the non-specialist 
reader. Let us define again, for the sake of completeness, the two relevant spaces: 1Z = M D 
and the order parameter space O as the space of real unit vectors v = v(xi, X 2 , ■ ■ ., xp) in 
the D dimensional Euclidean space E D with components denoted v r At the non-singular 
points of v, we also define the differential D — 1 form 

uj = Vj x dvj 2 A dvj 3 A • • • A dvj D e^ 2 '"^ D (Al) 


where e is the rank D totally anti-symmetric tensor. Finally, we define the topological index 
of the volume V as the quantity 


r(D/2) r 

" = WfA 

where it is assumed that uj is defined everywhere on the surface S surrounding the volume 
V. The normalisation comes from the formula for the surface area of ) the unit sphere 
in D — 1 dimensions, with T(x) being the gamma function. The quantity // is a topological 
invariant of the volume enclosed by surface S. 

Because of the normalisation condition, d{vjVj) = 2 Vjdvj = 0, at any non-singular point 
there is some non zero vk such that dry = A- ^2j^ k v jdvj. Hence at non-singular points, 


duj = dvj x A dvj 2 A dvj 3 A • • • A dvjj^e^ 2 "'^ 0 
= D\dVj 1 A dvj 2 A dvj 3 A • • • A dVj D 

= 0 . 


The first line follows because dd = 0, the second from the antisymmetry of the wedge 
product, and the third from the aforementioned fact that one of the components dv k is not 
independent and dvj Advj = 0. If the region V contains no singular points, then uj is defined 
everywhere in V and from the generalised Stoke’s theorem 


However, when the region does contain one or more singular points, the form c o (from Eq. [Tj) 
is not defined everywhere in the region and we must be more careful. Let us focus on some 
singular point p and denote its neighbourhood V p and S p = dV p , where we assume that the 


16 





surface S p contains no singular points. The space of unit length vectors in D dimensions is 
isomorphic to the D — 1 dimensional unit length sphere so the vector field v defines 

a mapping: 

S p s {D - 1] 

Because by assumption S p contains no singular points, it is homotopic to S (D ^ 1> , hence the 
above mapping is classified by the homotopy group = Z. The meaning of 

this is that the quantity p takes integer values and it is invariant to deformations in the 
region S p that do not cross singularities; i.e. it is a topological invariant. 

If we choose the neighbourhood V p small enough then the only singularity will be p and 
the index characterising it is 

^ = | p weZ -{°}' 

Appendix B: Proof for the position of the centroids when the system has N-fold 
rotational symmetry 

If the weights are solely a function of the intensity value, then an intensity pattern with 
rotational symmetry (cylindrical or discrete rotational symmetry) will yield a point at the 
symmetry axis. Defining the symmetry axis as the origin, this means p\ = 0. We will show 
this for N-fold rotational symmetry. In order to simplify the notation we will drop the 
dependence of the intensity I(yi, 1 / 2 ', Xi, x 2 ) on the position of the object (xi,x 2 ) and will 
drop the subindex of the vector p t , as the proof is the same for all vectors. 


P 



h (/ (y lf y 2 )) ( 2 / 12/1 + 2 / 22 / 2 ) d//i dy 2 

2 lT POO 


h (/ (r cos ( 6 ), r sin (0))) (r cos ( 9) yi + r sin ( 9) y 2 ) rdr d 6 


j o Jo 
N-l 


E 


n=0 1/0 J0 


h / r cos ( 9 + 


2nn\ 


, r sin ( 9 + 


27m 


x ( rcos [e + ^] yi + r sin ( 9 H—— 1 y 2 ) rdr d 9 


N~ 1 .27T 


-E 

n =0 


NJ \ N 

h (/ (r cos (9 ), r sin (9))) 


/o Jo 
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(Bl) 


X 


cos 



yi + r sin 




rdr d 9 


= 0 . 

In the second equality we have performed a change of coordinates: (yi, y 2 ) = (r cos (9),r sin(d)) 
In the third equality we have divided the azimuthal integral in N different parts. In the 
fourth equality we have used the assumed N —fold rotational symmetry of the intensity 
distribution. The last equality is just a result of the summation over the sin and cos 
functions. 
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